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Abstract 
We  investigate  the  complexity  of  the  problem  of  coordinating  the 
motion  of  many  planar  rigid  objects  in  a  two  dimensional  space.  Each 
object  has  a  simple  geometrical  description  (it  is  a  box  or  a  circle). 
The  number  of  objects  to  be  moved  is  part  of  the  problem  size.  The 
plane  where  objects  move  is  partitioned  into  two  regions:  The  region  of 
forbidden  points,  being  a  finite  union  of  areas  each  bounded  by  a 
polygonal  line,  and  its  complement.  During  motion,  objects  are  not 
allowed  to  overlap  or  to  intersect  the  region  of  forbidden  points,  but 
they  are  allowed  to  be  in  contact  with  the  boundary  of  obstacles  or 
other  objects.  We  analyze  several  variants  of  the  above  problem  of 
motion  planning  and  prove  that  they  are  NP-hard,  in  the  strong  sense. 
We  also  consider  simple  combinatorial  planar  graph-and-pebble  models 
which  capture  certain  non-geometric  aspects  of  coordinated  motion.  We 
provide  an  algorithm  for  planning  the  motion  of  k.  pebbles  in  an 
n-vertex  planar  graph,  which  is  0(n  ).  We  also  show  that  coordination 
of  motion  of  n-1  distinct  pebbles  in  an  n-vertex  graph,  from  a  set  of 
n-I  distinct  initial  to  n-1  distinct  final  positions,  can  be  decided 
and  planned  in  polynomial  time  in  n.  Finally,  we  show  that  certain 
problems  of  planning  the  motion  of  many  pebbles,  when  the  motion  is 
nonsymmetric,  or  when  the  final  positions  are  required  to  satisfy  a 
combinatorial  property,  are  NP-hard  in  the  strong  sense. 
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1.  Introduction 
1.1.  The  problem 

This  paper  examines  the  complexity  of  the  problem  of  planning  the 
motion  of  many  simple  rigid  objects  in  the  plane.  The  objects  are 
either  boxes  or  circles.  The  number  of  objects  to  be  moved  is  part  of 
the  problem  size.  The  plane  where  objects  move  is  partitioned  into  two 
regions  :  (a)  The  region  F  of  forbidden  polncs,  being  a  finite  union  of 
(not  necessarily  convex  or  finite)  areas  F,  ,  •eacVi' ot  which  is  bounded 
by  a  polygonal  line.  (b)  The  compleraeni:  of  F,^^  called'  thu  region  of 
legal  points.  A  position,  p(0^),  of  an  object  0  •  is  a  valid  position 
if  the  area  of  p(G  .)  n  F  is  zero,  and  also  if  the  area  of  p(O^)  ^  p(O^) 
is  zero  for  each  k  *  j  and  where  pO^)  is  the  position  of  object  9^  at 
some  instant  of  the  motion.  (i.e.  we  allow  objects  to  come  in  contact 
with  walls  -  the  boundaries  of  F  -  or  other  objects  during  motion,  but 
we  do  not  allow  them  to  be  able  to  enter  into  a  wall  or  overlap  with 
other  objects).  In  some  of  the  problems  to  be  examined  F  can  be  empty 
or  nnay  have  a  particularly  simple  geometric  description  (e.g.  the 
number  of  the  polygonal  lines  is  a  small  constant).  In  addition,  we 
are  given  an  initial  set  I  of  valid  positions,  I  =  {P^:  P^  is  the 
initial  position  of  object  Q ^}  and  a  set  of  sets  of  final  positions  FP 
=  {S.:  0.  is  an  object  and  S.  is  a  set  of  valid  positions  for  object 
0  .} .  In  some  of  the  problems,  S.  consists  of  finite  collection  of 
positions;  in  other  problems,  S.  is  defined  as  the  set  of  points 
surrounded  by  some  polygonal  line.  Another  way  to  specify  the  set  FP 
would  be  to  specify  US.  as  the  set  of  points  surrounded  by  a 
specific   polygonal  line  and  to  require  that  the  final  positions  of  the 
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objects  have  some  additional  property,   e.g.    cover  a  given  set   of 
points. 

The  general  two-dimensional  motion  planning  problem  is  the 
following: 

(a)  Given  (F,0,I,FP)  where  0  =  {0^,...,0^},  I  =  {Pp...,P^}, 
FP  =  {Sp  ...  ,S^}.  and  F  =  {Fp...  .F^^}  can  0  be  moved  from  I 

to  FP  via  a  c^ntl^uoy^.  mol ion  of  objects  so  that  ail  intermediate 
positions  ar^(.v^l;L4?  r   ^ 

(b)  Furthe^Taore-.,j_fLfg,the  •answer  to  (a)  is  yes,  is  it  possible  to  specify 
the  motion?    -  ^,  :  , 

The  input  length  of  an  Instance  J  of  the  problem  (a)   and   (b)   is 
defined   to  be   the  number  of  symbols  in  the  description  of  J  obtained 
from  some  fixed  encoding  scheme  for  (a)  or  (b),  applied  to  J. 
1.2.   Previous  work. 

Reif  [R-79]  examined  the  problem  of  moving  hinged  objects  in 
3-dimensionai  space  in  the  presence  of  obstacles.  In  this  work,  the 
shape  of  the  object  and  of  the  obstacles  was  allowed  to  be  complicated. 
The  problem  was  proved  to  be  P-SPACE  hard.  Hopcroft,  Joseph  and 
Whitesides  [HJW,82]  examine  the  problem  of  planning  the  motion  of  a 
complicated  linkage  in  the  plane,  and  show  that  this  problem  is  P-SPACE 
hard,  even  in  the  absence  of  obstacles.  Recently,  Hopcroft,  Schwartz 
and  Sharir  [HSS,83]  proved  that  the  coordination  of  motion  of 
complicated  "packed"  non-convex  bodies  in  a  2-dimensional  convex  space 
is  P-SPACE  hard.  They  remark  that  the  result  could  be  extended  to 
complicated  convex  bodies.  No  previous  research  attacked  the  problem 
of  coordination  of  motion  of  very  simple  objects  in  2-dimensions. 
(E.g.   disks,  boxes). 
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1.3.   A  combinatorial  view  of  the  problem:  Graphs  and  pebbles. 

Motion  planning  for  many  independent,  simple  objects  has  to  take 
into  account  (a)  The  geometry  of  the  obstacles  and  objects  and  (b)  The 
relative  positions  of  objects  and  obstacles,  which  forces  some  degree 
of  coordination  of  the  object  movements.  We  here  introduce  a 
combinatorial  model  which  captures  only  the^'^  ^problem  of  motion 
coordination,  stripped  of  its  geometric  aspects^.  -■' 

Consider  a  planar  graph  (or  digraph)  G(D)  'of  vertex  set  V  and  edge 
set  E  (arc  set  A).  Let  |v|  =  n  and  let  IC  V  Be  a  set  of  vertices, 
such  that  i  e  I  is  occupied  by  a  pebble  p^  (i.e.  an  1-1  assignment  of 
names  to  vertices  of  I).  In  addition,  a  collection  S  =  {Sj^ , . . .  ,S  i  ■[•  i  } 
is  given,  of  sets  of  vertices  of  V.  A  pebble  p,  currently  located  at 
vertex  i,  can  move  to  vertex  j  iff 

(1)  {i,j}  e  E  (or  (i,j)  e  A  in  case  of  a  digraph)  and  (2)  no 
pebble  is  located  at  j  before  pebble  p  moves.  We  call  such  a  move  a 
legal  move. 

The  combinatorial  motion  coordination  problem  is  as  follows: 

Given  (G  or  D,  I,  S),  is  there  a  way  of  moving  one  pebble  at  a 
time  (by  a  legal  move)  so  that  the  final  position  of  pebble  p^ 
(originally  placed  at  i)  belongs  to  S^  (for  i  =  1 ,2, . . . , 1 1 1 )? 

In  this  set-up,  unidirectional  motions  are  modelled  by  directed 
edges  and  we  also  may  consider  one-pass  vertices :  A  vertex  is  one-pass 
if  no  additional  pebbles  can  visit  it  as  soon  as  a  single  pebble  passes 
through  it.  In  this  paper,  we  provide  evidence  that  the  (geometric) 
motion  planning  problem  is  harder  than  the  problem  of  motion 
coordination  of  pebbles  in  planar  graphs.   We  give  an  example  in  which 
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the  two  problems  are  equivalent:  (Fig.  1 )  In  part  A,  the  set  of  legal 
points  is  an  orthogonal  ABCD  which  initially  has  nm-1  squares  (objects) 
in  it,  and  a  free  space  which  can  fit  exactly  one  object.  In  part  B,  a 
lattice  graph  of  nm  vertices  contains  nm-1  pebbles  and  a  free 
(unoccupied)  vertex.  We  can  define  an  1-1  correspondence  between 
objects  and  pebbles,  and  between  initial  and  final  positions  of  objects 
and  pebbles.  Any  ^a^qu^jcui^  of  valid  moves  for  one  of  the  problems  A,B 
corresponds  to  a  valid  -fOCjJtic^n  planning  for  the  other  case.  Problem  A 
is  a  genera;iza,ti,C5iv,,of  ^tli« -.ISTpuzzle  (see  [S,  1879]).  Problem  B  may  be 
generalized  to  a  puzzle  based  on  an  arbitrary  graph.  In  this  paper  we 
show  that  the  decision  problem  and  the  motion  planning  problem  for  the 
genrealized  puzzle  B  (and  hence  for  A)  can  be  solved  in  polynomial  time 
in  the  number  of  .vertices  (squares  for  part  A). 
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2.  NP-hardness  results  for  geometric  motion  planning  problems   in  the 
plane. 

This   paper  shows   the   following  problems   to  be  NP-hard  in  the 
strong  sense: 

BOX  MOTION  PIANNING 
(PI)  A  collection  of  n  boxes  is  given,  each  specified  by  its  two  sizes 
(xj^.y^),  i  =  l,2,...,n.  A  set  of  initial  positions  I  =  {ip...,i^}  is 
given,  where  each  i,  specifies  the  center  of'  b?)'x'0^  and  the  direction 
of  one  of  its  edges.  A  set  of  final  positioWs"F?  ==  {lp...,f^}  is 
given,  one  for  each  object,  where  each  f .  specifies  the  final  position 
of  the  center  of  G  .  and  the  final  direction  cf  one  of  its  edges. 
Question:  Can  we  move  the  boxes  from  I  to  FPV 

Remark:  Reduction  from  3-PARTITION.  Objecc  sizes  can  be  restricted  to 
be  at  most  0(n").  We  allow  the  final  positions  of  all  but  one  of  the 
boxes  to  coincide  with  the  initial  positions.  Object  sizes,  coded  in 
binary,  are  part  of  the  input.  The  problem  remains  rJP-hard  (in  the 
weak  sense  however)  even  if  the  number  of  lines  in  polygons  is  a  small 
constant.   In  that  case,  a  reduction  from  SUBSET  SUM  can  be  used. 

BOX  MOTION  PLANNING  WITH  'DON'T-CARES' 
(P2)  A  collection  of  n  boxes,  obstacles  and  initial  positions  I  is 
given,  as  in  (PI).  A  set  FP  =  {S,,...,S  }  of  sets  of  final  positions 
(one  set  for  each  object)  is  given.  Some  of  the  Si's  may  be  singletons 
(uniquely  specified  positions)  while  some  others  are  areas  whose  union 
is  specified  by  polygonal  boundaries. 

Remark:  Transformation  from  3-PARTITION.  The  same  remarks  as  in  (Pi) 
apply.   The  numbr   of   objects  whose   final  position   is  completely 
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speclfled  may  be  any  small  constant  >    1.  The  polygons  can  be  restricted 
to  be  a  collection  of  horizontal  and  vertical  line  segments. 
(P3)  DISKS  MOTION  PLANNING 

A   collection   of   n   disks   of   various  radii,  r . ,  i=l,2,...,n  is 
given,  together  with  a  set   of   polygonal   obstacles.    A  set,   I,   of 
initial   positions,  one  fixed  position  for  each  disk,  is  given.   A  set, 
FP,  of  final  posltiprii^,  op^^  fixed  position  for  each  disk,  is  given. 
QUESTION:  Can  we,_  move- tji^,  disk^,  from  I  to  FP? 

Remark:  Transf  oxmatipn^  frogi.,  3--PARTITI0N.   We  allow  the  final  positions 
of  all  but  one  of  tKe  disks  to  coincide  with  the  initial  positions. 
(P4)  DISKS  MOTION  PLANNING  WITH  'DON^T-CARES^ 

A  collection  of  n  disks,  obstacles  and  initial  positions  is  given, 
as  in  (P3).  A  set  FP  =  {Sp...,S^}  of  sets  of  final  positions  (one  set 
for  each  object)  is  given.  Some  of  the  S.'s  may  be  singletons 
(uniquely  specif  led.  positions )  while  some  others  are  areas  whose  union 
in  specified  by  polygonal  boundaries. 
QUESTION:  Can  we  move  the  disks  from  I  to  FP? 

Remark:   Transformation  from  3-PARTITION.   The  number  of  disks  whose 
final  position  is  completely  specified  can  be  any  small  constant  >  1. 

DISK  MOTION  COVERING 
(P5)  Given  n  disks  of  the  same  radius  (r)  and  n   initial  positions   of 
their   centers.    Can  we   move   them  in  the  plane  so  that  their  final 
positions  cover  a  given  set  of  m  points? 

Remarks:  Remains  NP-hard  even  if  r  is  constant  and  no  osbstacles  exist. 
Reduction  from  3-SATISFIABILITY 

Note:  References  for  the  NP-coraplete  problems  which  furnish  a  basis  for 
our  arguments,  can  be  found  in  [GJ,  79]. 


3.  NP-hardness  results  for  planning  pebble  motion  on  a  graph 

ONE  PASS  PEBBLE  MOTION  IN  PLANAR  GRAPHS 
(P6)  Given  a  planar  graph  G  =  (V,E)  with  ail  vertices  in  V  being 
one-pass.  Given  also  a  collection  I  =  {s,,...,s.  }  E  V,  s  = 
{tp...,t^}  E  V  and  ^i.'^j  such  that  j  *  i  L^  *  l^,  s^  i^  Sj,  s^  #  ty  A 
pebble  named  p.  is  originally  placed  on  each  s^,  i  ~  1 , . ^ . ,k.  Is  there 
a  way  of  moving  each  p.  from  initial  pcr.ltlon's,  to  frlnrl  position  t. 
by  a  sequence  of  legal  moves?  Re ma rk  Transf pncation  from  "k-vertex 
DISJOINT  PATHS  IN  PLANAR  GRAPHS",   k  is  p^.t  of  the  input. 

PEBBLE  MOTION  WITH  NON-ADJACENCY 
(P7)  Let  there  be  given  a  planar  graph  G  =  (V,E)  such  that  maximum 
degree  of  each  vertex  u  e  V  is  <  4,  together  with  an  initial  position 
I  =  ^1 .  •  •  •  >ii<;}  5.  V  of  k  pebbles.  Given  a  set  of  final  positions 
FPC  V.  Can  we  coordinate  the  motion  of  pebbles  so  that  the  final 
positions  have  the  property  that  no  two  pebbles  are  neighbors  and  ail 
final  positions  belong  to  FP? 
Remark   Reduction  from  INDEPENDENT  SET  IN  CUBIC  PLANAR  GRAPHS. 
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4.      Proofs    of    the   NP-hardness    results 
4,1    Proof   of    (PI) 

We  shall  transform  the  3-PARTlTION  problem  to  the  BOK  MOTION 
PLANNING  problem.  The  3-PARTITION  problem  [G7,  pg.  96]  is  defined  as 
follows: 

A  finite  set  A  is  given  of  3m  elements,  a  bound  B  e  Z"^,  and  a  size 

s(a)  e  Z"*"  for  each  a  e  A,  such  that  s(a)satisf ies  B/4  <  s(a)  <  B/2   and 

such   that  1    is(.a)  i  =  -raBc'   Can  A  be  partitioned  into  m  disjoint  sets 

acA 
S,,S2,...,S   such  tliaC.for-^l  <      I     <.      m,      I        s(a)=B?   Given  any 

ae  S  J 

instance  I  of  the  above  problem,  we  shall  provide  a  polynomial  time 
transformation  to  an  Instance  F^  of  the  BOX  MOTION  PLANNING  problem,  as 
follows:  (Fig.  2).   Let  n  =  3m  . 

The  elements  a^  of  A  will  correspond  to  small  boxes  a^,...,a^  each 
rt: 

placements  of  all  a.  are  within  the  orthogonal  ABCD  of  width  CD  =  Xq  > 
B  and  height  BC  >  2m+l.  To  the  left  of  ine  AD,  there  are  in  slots, 
each  of  which  has  length  ED  =  B  and  height  w  =  1.  Successive  slots  are 
separated  by  obstacle  walls  of  height  equal  to  1.  The  object  0,  must  be 
moved  so  that  its  center  goes  from  I,  to  F.  Its  length,  i,  is  at  least 
2AB.  Its  width  is  equal  to  1  (equal  to  MN  or  PQ  in  Fig.  2).  The 
vertical  size  of  O2  is  at  least  BC.  In  order  for  0-^  to  pass  through, 
O2  must  come  down,  in  complete  contact  with  the  walls  BC,  CD  and  the 
line  AD.  Finally,  the  space  ABCD  is  such  that  an  arrangement  of  all 
objects  a.  in  it  is  possible,  as  in  Fig.  3.  This  arrangement  allows 
objects  a,,..., a  to  "rotate"  into  the  space  of  ABCD  so  that  any  object 
a.  can  be  at  any  of  the  slot  entries.  This  allows  any  possible 
arrangement  of  objects  into  the  slots  (given  the  size  constraints). 
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Any  motion  planning  which  succeeds  to  move  Oj^  from  I^^  to  Fj^  must 
at  some  point  in  time  have  0^  in  contact  with  CD.  This  is  only 
possible  if  the  boxes  o  ■,,...  yO  can  be  "packed"  into  the  m  slots.  Note 
that  each  slot  can  fit  at  most  3  boxes  since  each  sCa^^)  >  B/4.  Also, 
if  one  slot  gets  only  two  boxes,  this  will  imply  that  another  slot  will 
get  four  boxes,  which  is  impossible.  So,  each  slot  has  to  get  exacty 
three  boxes.  Now,  if  there  is  at  least  '  one"  s'lot  which  is  only 
partially  filled  (i.e.  the  sum  of  the  horizontal  sizes  of  the  3  boxes 
into  the  slot  is  less  than  B),  then  there  must"  be 'another  slot  whose 
boxes  have  sum  of  horizontal  sizes  greater  than  B.  (This  is  so,  since 
\  s(a,)  =  mB).  The  above  conclusion  cannot  be  realized,  since  slots 
are  of  length  B  and  object  O2  blocks  the  entrance  of  the  slots.  Hence, 
object  0^  can  pass  through,  if  and  only  if  the  assignment  of  objects  c^ 
to  slots  is  such  that  their  horizontal  sizes  a(aj.)  solve  the 
3-PARTITION  problem. 


the  above  construction  works  even  if  the  maximum  object   size   (AN)   is 
0(n5). 

Note  that  a  similar  transformation  can  be  done,  from  SUBSET  SUM, 
by  using  only  2  slots.  This  shows  that  even  when  the  number  of  lines 
in  the  polygonal  boundraies  is  constant,  the  problem  is  still  NP-hard 
(in  the  weak  sense  however),  due  to  very  large  object  sizes  and  the 
economic  (binary)  integer  representation. 

4.2.   Proof  of  (P2) 

We  will  use  a  transformation  from  3-PARTITION.   The  construction 
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ls much  like  Figure  2  (see  Fig.  4).  The  elements  of  A  will  correspond 
to  boxes  ^i.''».o  (n  =  3m),  of  height  1  and  length  s(a^).  Again,  the 
initial  position  of  object  O2  allows  all  a.'s  to  be  positioned  in  the 
space  ABCD  as  in  Fig.  3.  This  enables  any  motion  planning  algorithm  to 
form  any  assignment  of  '^  ]> ' ' '  '^n  ^^  the  m  slots  (given  the  size 
constraints).  The  final  position  of  O2  is  in  complete  contact  with 
walls  CD,  BC  and  line  AD.  The  final  positions  of  a^,...,a^  all  belong 
to  the  area  of  the  m  slots.  All  objects  reach  their  final  positions 
iff  the  s(a.)'s  solve,  the  3-PARTITION  problem  (similar  arguments  as   in 

• 
Remark.   The  motions  allowed  in  (PI)  and  (P2)  can  be  restricted  just  to 
vertical  and  horizontal  motions.   The  above  results  hold  also  for  Oj^'s 
being   ellipseso     Similar   remarks   as   in   4.1   about   a  possible 
transformation  from  SUBSET  SUM. 

4.4.   Proof  of  (P3) 

We  will  use  a  transformation  from  3-PARTITION.  Given  any  instance 
I  of  the  3-PARTITION  problem,  we  will  transform  it  to  an  instance  F-j.  of 
DISKS  MOTION  PLANNING,  in  polynomial  time,  so  that  the  answer  to  I  is 
YES  iff  the  answer  to  F^  is  yes.  Let  n  =  3ra.  Cosider  the 
configuration  of  Figure  5. 

Our  configuration  consists  of  two  large   circles,   C^,   Cy  and  n 

small   circles  a  ^, . . .  ,a ^.      The  radii  of  the  circles  are  r^^,  r  ,  r^  . . .  ,  r^ 

s(ai) 
respectively,  with  r.  =  1  +  c-,  where  e^^  is  exactly  — ^5— 7-*  where  S  =  mB 

and  a.  t    A.  It  is  clear  that  the  representation  of  each  e  ^^  is  linear  to 

the  representations  of  s(a.),  n  and  B.  The  boundary  G   is   part   of   a 
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canonical  polygon,  of  sufficiently  large  number  of  vertices,  so  that  it 

can  be  superscribed  to  C^  and  at  the  same  time,  when  Cy  is   in  contact 

with   G,   it  is  impossible  for  any  of  the  objects  o ^,    i   =  l,...,n  to  be 

anywhere  but  in  one  of  the  m  "slots"  Tj^,...,T^.   (See  Figure  6.)  For  G 

to  have   such  a  property,  it  is  enough  to  select  c^   so  that  (1)  y  >  m 

and  (2)  E'Y'  -  EY'  «  2  (Fig.  6).   Since  E'Y'  -  EY'  =  i?!J_  and  EE '   = 

^Y 
ry  "^Y      9  -1 

—  •    TT  ,   we  have  E'Y'    -  EY'    =  — it  ^ ,    so   it    is    enough   to   sslect    ry  =   kir 
^  k^  , 

and  k  >    n^S^,    in  which    case   EE'    =    1    and  E'Y'    --E^'-jA  P^-^]    «   2      (Eq. 

n  S 
*).      Furthermore,    (1)    is   satisfied,    since  k  >,  n  s4,J>   ^r.  =   2in  for  n  >    2. 

Let  e^j^^  =  max{e^},  i  =  l,2,...,n.   The  distance  w  (width  of  each 

slot)  is  2  +  2£^^^   and  the  distance  I    (length  of  each  slot)  is  £  =  6  + 

1   "  '         1 

-  I      e,.    Since  e,      <  — ,   it   is   clear   that  at  most  three  of  the 

^   1=1  n2 

a, a   can  fit  in  each  slot,  if  their  centers  lie  en  the  same   line. 

1 '    '  n 

The  same  is  true,  even  if  one  tries  to  "pack"  disks  into  the  slots  as 
much  as  possible.  The  reason  is  that  the  gains  in  occupied  space,  when 
the  disks  are  not  concentric,  are  quadratic  on  the  Ej^'s  and  hence  the 
sum  of  any  three  of  such  gains  is  much  less  than  min{e^}  (see  Fig.  7). 
From  Figure  7,  the  gain  in  space  by  packing  a^  ,  a ^,  so  that  their 
centers  are  not  on  a  line  parallel  to  the  slot's  direction,  is 


(LL')2  .  C2e^,-e,-C2)\^^^.   ^ 


KL  +  KL'  2  t°ax 

Hence,  the  total   gain  by  packing   three   disks   as   economically   as 
possible  is 
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_       ^   max~s(aj)        , 
0(3  ^2^,)  =  0(3     ^/  ]  <  O(^)  <  e^^  .  (Eq. 

n  S  n  S 


In  order   for  C^  to  move  from  X  to  X',  Cy  has  to  move  down  to  come  in 

contact  with  G.  At  this  moment,  the  objects  a^,...,a^  have  to  split   in 

m  sets   of  cardinality  three  each  and  move  into  the  ra  slots.   (If  less 

than  three  of  th-e  a  . '^ -mbVe  "into  some  slot,  then  more  than   three  will 

have  to  move  "into  anotHeir"'slot ,  but  this  is  impossible.) 

This   implied'  tha'e"  the  re  exist  S^,...,Sj^,  S^C  A,  such  that  \5^\    = 

3,  for  i  =  l,2,...,m  and  Vi,  if  S.  =  {a,  ,a,.  ,a^  }  then  a,-  ,a,-  ,0^   are 

1     1^'  i2'  13        ii'  I2'  I3 

in  the  slot  T^.    -  ' 

We   now  claim  that  ¥i   =   I,2,...,m,  T    sCa^)  =  B.  To  show  it, 

a^eSi 

suppose  (for  contradiction)  that  there  exists  j:  I  ^^^k^  ^     ^*   ''^hen 

,   V     ^k=  ^  i 
it  must  be  the  case  that  there  exists  j':  )    s(av.)  >  B  (by  at  least  1 

unit).   Let  A  =  I  s(a,^)  -  B.  Then  A  >  1  and  6  =  -A.-  >  -^tT'        '^^^ 

a^eS.  n-S^     n^S^ 

three  disks  a^   (such  that  a^  e  S.)  have  to  be  "packed"  into  the  slot  T. 

in  such  a  way,  to  gain  at  least  a  5  length  (because,  if  we  pack  them  so 

that   their  centers  are  on  a  line  parallel  to  the  direction  of  T.,  then 

they  will  not  fit  into  T.,  by  Eq.   *).   But  we  showed   that   the   total 

gain   by  packing   3  disks  as  economically  as  possible  is,  by  Eq.   **, 

less  than  0(_ )  which  is  <  5 .   This  implies  that  the  three   disks   of 

T.  will  not  fit  into  the  slot,  and  that  contradicts  our  assumption  that 

Cy  is  in  contact  with  G.  We  conclude  that  Vi  =  l,...,m,  I  s(a^)  =  3. 

Conversely,   if   this  last  equality  holds,  then  the  disks  o  -^y . . .  ,a^   can 

be  separated  into  ra  sets  of  three  disks  each,  such  that  the  sum  of   the 

diameters  in  each  set  is  6  +  =  6  + =  6  +  —  )   e  ,-  =  2. .   This 

„2c:2       m  „2<;2       m.  ,  ^ 
n  S  n  b         1=1 
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means  that  it  is  possible  for  the  disks  to  be  packed  into  the  ra  slots, 
and  to  allow  C^  to  be  in  contact  with  G  and  C^  to  pass  from  X  to  X'. 

e 
Remarks 

(1)  The   final  positions  for  C^,  a^,...,a^  may  just  coincide  with  the 

initial  positions.  "    ,■  ' 

s(a^) 

(2)  Any  choice  of  e  ^)'s   smaller  than       -wcrifM  do,  "'a:s   soon  as   their 

n  "S 
description   is  polynomially  related  to  that  c>f  sCa.) 's.   By  our  remark 


would  use  nonstandard  analysis  and  "ideal"  infinitesimals,  used  as 
indetermi nates  in  some  of  the  calcuations  (see  [Davis,  77]).  The 
"ideal"  infinitesimals  then  solve  the  3-PARTITION  problem. 

4.4.  Proof  of  (P4) 

In  (P3) 's  construction,  allow  the  objects  Cy,  ai,...,a  to  be 
'care-frees'  (i.e.  their  final  positions  are  only  restricted  to  be 
legal).   Use  the  same  proof. 

4.5.  Proof  of  (P5) 

Our  construction  is  similar  in  flavor  to  [FPT,  81].  See  Appendix 
I  for  a  detailed  argument. 


4.6.  Proofs  of  (P6),  (P7) 

See  Appendices  II,  III. 
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5.   An  0(n  )  algorithm  for  motion  planning  of  k  pebbles  in  an  n-vertex 

planar  graph. 

Let   (G,I,S)   be  an   instance   of    the   combinatorial   motion 

coordination  problem  with  G  =  (V,E),   |V|   =  n,  |l|  =  |S|  =  k.  (Same 

algorithm  in  case  of  directed  graphs).    Let   us   define   a   new  graph 

G   =  (V  ,E  )   where   V  =  V^   (i.e.   VxV.  ..xV,  k  times)  -  X,  where  X  is 
m     m'  m^  m 

the  set  of  k-tuples  .  (y,  , . . .  ,vi^) ,  v.  e  V,  i  =l,...,k,  and  3  i,j,L*2 
such  that  V,  =  v-.  Two  vertices  x,y  of  V  are  connected  by  an  edge 
eG  E^  iff  X  =  (u^,...»,U|^)  y  =  (v^,...,Vj^)  and 

(1)  3-i,  i  <  i  <  k  such  that  {u^^.v^}   E  and 

(2)  u.  =  V.   for   j  *  i  and 

(3)  Uj  7t  v^  for   j  *  i 

Intuitively,  x  represents  the  position  of  the  k  pebbles  before  the 
move  of  one  pebble  and  y  the  position  after  the  move.   It  is  clear  that 

Remark:  Any  motion  coordination  algorithm  defines  a  path  in  G^, 
connecting  I  C  V^  to  S  C  V^^^. 

Hence,  we  can  plan  the  motion  of  pebbles  from  I  to  S  by 
constructing  G  and  testing  if  I  and  S  are  connected  in  G^.  These 
operations  have  worst  case  time  0(|E  |)  and  the  space  needed  is  again 
0(|E^|)). 

We  now  will  state  a  bound  on  |E  |  based  on  the  planarity  of  G. 
(The   obvious   bound   is   0(|V^|2)   =  0(n'^'*^)).    For   any   x  £  V^,    if 

x  =  (uj^, . . .  ,uj^)  then  the  degree  of  x  in  G^^,  call  it  d^ix)    is  bounded  by 

k 

II      d(uj)  where   d(u^)    =  degree    of   u^    in  G.    But 
i=l 

k 

K\  =y    I      V^)  =1    I        n    d(ui) 

xeV  xeV„   i  =  l 
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Since   G   is   planar,   we  have 


^      d(u)    =  0(n)    implying      [    I      d(u))^  =  OCn'^)    . 


^■-^  at.V 


But,    since   d(u)   >    0  ¥u  V  we   have 
k 


I  J[      d(ui)  <  (    I      d(u)]^   =  OCn^^)    ..,^,^.    „_   -,^ 

xtV^  1=1  utV 

Hence  •' -■   g    ,-    iEn.l.lff' 

|E J  =0(n^)    .  ■"''•    ""^       ''■' 
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6.  Planning   the   motion   of   n   pebbles   in   an   n+1-vertex   graph   is 
polynomial  time. 

Let  G  be  any  finite  simple  graph  with  vertex  set  V  of  cardinality 
n+1.  By  a  "pebbling"  we  mean  the  placement  of  pebbles  l,2,...,n  on 
distinct  vertices  of  G,  leaving  one  vertex  unoccupied.  Formally,  a 
pebbling  on  G  is  the  bijective  mapping  f:  V  ->■  {  1 ,2, .  . .  ,n,0}  .  The 
vertex  x  with  xf  =  0  is  said  tc  be  unoccupied  in  f .  (Ail  mappings  will 
be  written  on  the  right  and  composition  of  mappings  gose  from  left  to 
right.)  Two  pebblings  f,  g  on  G  are  adjacent  iff  g  can  be  obtained  from 
f  by  moving  a  pebble  along  an  edge  of  G  onto  the  vertex  which  is 
unoccupied  in  f. 

The  relation  of  adjacency  defined  on  the  pebblings  is  symmetric 
and  irreflexive  and  so  it  defines  a  new  simple  graph,  G  with  vertex 
set  V  consisting  of  all  pebblings  of  G,  two  pebblings  being  joined  by 
an  edge  in  G^  iff  they  are  adjacent.   Let  the  graph  9 g  be  as  in  Figure 

7.  R.  Wilson,  [W,  74],  proved  that: 

Theorem  6.1  [W,  74],  Let  G  be  a  finite  simple  nonseparable  graph  other 
than  a  polygon  or  9 g.  Then  G  is  connected,  unless  G  is  bipartite,  in 
which  case  G  has  exactly  two  components.  In  this  latter  case, 
pebblings  f,  g  on  G  having  unoccupied  vertices  at  even  (respectively, 
odd)  distance  in  G  are  in  the  same  component  of  G  if  and  ony  if  fg~ 
is  an  even  (respectively,  odd)  permutation  of  V(G). 

Corollary  6.1.   The  decision  problem  for   the   coordinated  motion  of 
n-pebbles   in  an  n+1-vertex  graph  from  a  fixed  Initial  to  a  fixed  final 
position,  is  polynomial.   (In  fact,  it  can  be  decided  in  0(n'-)  time.) 
Proof:  Just  test  if  G  is  bipartite.   If  not,  then  answer  'yes'.    If   G 
is   indeed  bipartite,  then  find  the  distance  of  the  unoccupied  vertices 
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in  the  Initial  (f)  and  final  (g)  positions  of  the  pebbles  and  check  for 


6.1. 


(The  case  of  9 q  can  be  tested  in  constant  time  and  the  case   of   a 


polygon  is   testable  in  linear  time,  since,  given  an  initial  pebbling, 

I,  of  the  polygon,  the  only  reachable  pebblings  are  the  pebblings  which 

one   can  get   from  I   by  moving  each  pebble   of  I  clockwise   or 

counterclockwise  through  legal  moves.)  c 

Corollary  6.2.   Let  G  be  not  bipartite,  polygon  or'  the  graph  9q.   Then, 

it  is  possible  to  move  k  pebbles  from  any  inlti'ar  to  any  final  position 

(I  <  k  <  |V|  -  1). 

Proof:  By  Theorem  6.1,  if  we  extend  the  initial  and  final  pebblings  (in 

any  arbitrary  way)   by  adding  new  pebbles  so  that  the  total  number  of 

pebbles  is  |v|-l,  then  the  graph  G^  is  connected.   This   implies   that 

the  graph  G^  ,  for  the  k  pebbles,  is  connected. 

We  now  proceed  to  show  that  the  problem  of  planning  the  motion 
(i.e.  finding  a  squence  of  legal  moves)  for  n  pebbles  in  an  n+1  vertex 
graph,  is  polynomial  time. 

Driscoli  and  Furst  recently  proved  that  « 

Theorem  6.2  [DF,  83].  Any  group  represented  by  generators  that  are 
cycles  of  bounded  degree  has  O(m^)  diameter,  where  m  =  the  degree  of 
the  group.  I.e.  ,  the  longet  product  of  generators  required  to  reach 
any  permutation  in  the  group  is  O(m^).  Furthermore,  such  "short" 
products  can  be  found  in  polynomial  time. 

For  each   path  p   =   (xq,...,x^)   in  G   (|V(G)|  =  n+1)  define  a 
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Op  =  (XQXP  (X1X2)  (X2X3)  ...  (x^_ix^) 
(A  transposition  (xy)  is  a  permutation  on  V(G)  that  interchanges  x  and 

y.) 

Proposition  6 . 1  [W,  74].  Pebbiings  f,  g  of  G  are  in  the  same  component 
of  V  iff  f  =  a  g,  for  some  path  p  in  G  from  f~  (0)  to  g~  (0). 
Proposition  6.2  [W,  74]^.  ^  Fpr  each  vertex  x  of  G,  let  r(x,y)  be  the  set 
of  all  permutations  a  of  V(G)  where  p  is  a  path  from  x  to  y  in  G.  Let 
r(x)  =  r(x,x).  For  each  x,  T (x)  is  a  group  of  permutations  of  V(G), 
each  fixing  x.  If  p  is  a.  path  from  x  to  y  in  G,  then  r(x,y)  =  r(x)o  = 
OpPCy). 

Theorem  6 .3  [W,74].   Let  G  be  a  finite  simple  nonseparable  graph  other 
than  a  polygon  or  9 q.   Then  Vx  c   V(G), 

r(x)  =  sym(v(G)  -  {x}] 

(the  symmetric  group  of  all  permutations  of  vertices  V(G)  -  {x} )  unless 
G  is  bipartite,  in  which  case 

r(x)  =  ait(v(G)  -  {x}]  . 

(the  alternating  group  of  all  even   permutations   of   vertices  V(G) 

{x}).    Furthermore,   in   the  second  case,  the  generators  of  r(x)  are 

3-cycles. 

Theorem  6.4.   A  polynomial  in  length  sequence  of  legal  moves   for   the 

coordinated  motion  of  n  pebbles  in  an  n+1  graph  (other  than  a  polygon 

or  Gg)  can  be  found  in  polynomial  time  (in  n). 
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Proof :  Let  f,  g  be  Che  initial  and  final  pebblings  of  G  and  let  us 
assume  that  the  decision  problem  on  f,  g  has  the  answer  'yes'.  Let  p 
be  some  path  in  G  from  f'^cji)  to  g~^(.})).  Let  f  be  equal  to  o  g.  Both 
f,  f  are  permutations  of  V(G)  and  both  fix  the  same  "unoccupied" 
vertex  (the  f~  (<^)).  Hence,  f,  f  are  elements  of  the  group  r(f~  ((}))) 
(by  Prop.  6.2).  But  then,  by  Theorem  5.3,  r(f~^((j)))  is  generated  by 
bounded  degree  cycles  and  hence  (by  Theorem  6,2)  the  product  of 
generators  which  moves  from  f  to  f  can  be  foundf  in  poiynomial  time   in 
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7.  Discussion  and  further  work 

Any  pebble  motion  coordination  algorithm  has  to  find  a  path  from  I 
to  S  in  G^.  Hence,  the  diameter  of  the  worst  case  G^^^  is  a  lower  bound 
for  the  time  complexity  of  any  pebble  motion  planning  algorithm. 

Corollary  7. 1      ,  •.   !■  •  - 

If  the  planaif  grapJi.G,  in  an  instance  (G,I,S)  of  the  pebble  motion 
planning  problem,  , is  k-'connected,  then  motion  can  be  planned  from  I  to 
S  in  OCk^n^)  time. 

Proof.  Since  G  is  k-connected,  there  are  k  vertex-disjoint  paths  one 
for  each  pair  (s^.t^^),  s^  6  I,  t^  S  S.  The  algorithm  has  to  find  the 
paths  and  move  the  pebbles  along  them. 

We  remark  here  that  motion  planning  becomes  tricky  only  when  G's 
connectivity  is  less  than  k.  Further  work  will  have  to  do  with  the 
relation  of  the  complexities  of  the  geometric  and  combinatorial  motion 
coordination  problems.  We  believe  that  the  geometric  motion  planning 
problem  is  much  harder  than  NP,  even  for  simple  independent  objects 
like  disks  or  boxes  and  simple  obstacles  (orthogonals ).  The  problem  is 
conjectured  to  be  P-SPACE  hard. 

We  conclude  with  one  conjecture  and  one  open  problem. 
Conjecture  7. L.   There  is  a  planar  graph  G  =  (V,E)  and  positions   I,   S 
with   |l|   =  |S|  =  k,  such  that  the  number  of  legal  pebble  moves  needed 
to  move  ail  pebbles  from  I  to  S  is  at  least  n  ,  where  n  =  |v|. 
Open  problem  7. 1.   What  is  the  complexity  of  the   following:   Given  an 
instance   (G,I,S)   of   the  undirected  combinatorial  motion  planning 
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probiem,   with  G  =  (V,E),  G  planar  and  bipartite,  |v|  =  n,  |l|  =  |S|  = 
k,  k   varying  with  n,  can  we  move  k  pebbles  from  I  to  S? 
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Appendix  I 

Given  any  instance  I  of  2SAT  (M  variables,  N  clauses)  we  shall 
give  a  polynomial  time  translation  to  an  instance  F^  of  (P5).  The 
answer  to  (Pl)'s  instance  F^  will  be  yes  iff  the  3SAT  instance  I  is 
satisfiable.  Each  variable,  i,  corresponds  to  a  loop  of  points,  with 
the  following  property:  Each  pair  of  successive  poin*-salnng  the  loop 
can  be  covered  by  one  disk  and  any  given  di sic.  can  eovgr  jao  more  than 
two  points  along  the  loop.  This  is  so  e..7,.  if  the- distEnce  of  any  two 
successive  points  in  the  loop  is  d-e  ,  d  being  .thp-^di^k  diameter  and  e 
very  small  and  the  loop  is  as  in  Fig.   A.I,  with  AK.>  d- , 

Then,  each  loop  of  2k.  points  can  be  covered  bv  k.  disks,  in  one 
of  two  possible  ways,  one  of  which  will  correspond  to  TRUE  and  the 
other  to  FALSE  value  assigned  to  i.  The  initial  disk  position  with 
respect  to  each  loop  is  such  that  (1)  each  disk  can  move  to  an 
acceptable  final  position  without  being  prohibited  by  other  disks  and 
(2)  motion  of  disks  to  acceptable  final  positions  is  not  prohibited  by 
disks  already  placed  in  final  positions,  see  Fig.   A. 2. 

The  final  position  of  disk^  is  any  cardinality  2  subset  of  {X^, 
Af,  A^^.]^}  where  1,  i+1  are  taken  modulo  k^. 

Since  the  construction  must  be  embedded  in  the  plane,  we  have  to 
implement  crossovers.  At  crossovers,  the  loops  have  one  point  in 
common  and  they  can  share  one  disk.  When  sharing  one  disk,  the  value 
assignments  are  propagated  independently.   (Fig.   A.  3,  A.  4). 

For  each  clause  C.,  j=l,...,N  of  I,  there  is  a  region  in  which  the 
3  variable  loops  (which  correspond  to  the  3  literals)  come  to  close 
proximity  (let  us  call  them  L.^  L-^,  L.3).  In  the  clause  region  there 
are  4  additional  clause  points  (J,,  J2,  J3,  J^).   Each  J. (1=1,2,3)   can 
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be   covered   "for   free"  only  if  the  corresponding  literal  (to  the  loop 

^  Aj_)   is  TRUE.   Any  three  or  less  of  J^,    Jo,  J-j,  J4  can  be  covered  by   a 

clause   disk,   D.,   initially  placed   in   the   area   of  the  clause  and 

initially  covering  J^,  but  not  all   4   of   the   clause   points   can  be 

covered  by  the  clause  disk.   (Fig.   A.  5,  A.  6). 

Note  that'  the'  placement  (initial)  of  disk  D  does  not  prohibit  the 

motion  of  the''diskS^Of  ^ari^ble  loops.   The  final   positioning  of   the 

disks   of  vatiable-lob^s  mast  not  prohibit  the  motion  of  D    The  final 

position  of  D'.  is  ahy-sabset  X  of  {J^,  J2,  ^3,  J4}  with  I  <.     |X|  <  3. 
■'  N 

In  total,  we  hav^-'  ^   k^  +  M  disks,  all  of  the  same  diameter.   The 
1=1   N 
total  number  of  points  is'  2  J^   k.  -  N  +  4M.   All  points  can  be  covered 

1  =  1 
iff  at  least  one   literal  per   clause   is   true   (because,   then,   the 

corresponding   clause  point   is   covered   for  free  by  some  disk  of  the 

variable  loop  corresponding  to  that  literal).   QED 


Note:  We  took  care  in  the  construction  so  that  disks  can  be  moved,  one 
by  one,  from  their  initial  to  some  possible  final  position,  without 
their  motion  being  prohibited  by  the  initial  or  final  positioning  of 
other  disks,  in  any  other  way  than  restricting  the  set  of  final 
positions,  but  never  making  it  empty. 

No  obstacles  are  needed  for  the  transformation. 
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Appendix  II 


Proof  of  (P6) 

We  shall  transform  any  instance  I  of  k-VERTEX  DISJOINT  PATHS  IM 
PLANAR  GRAPHS  to  an  instance  F^  of  P6,  [GJ,  pg.  217],  so  that  P6  will 
have  a  YES  answer  iff  I  has  k  disjoint  paths.  Given  I,  we  just  put 
pebbles  in  the  vertices  s,.  The  fact  that  vertices  of  V  are  one-pass 
guarrantees  that  any  motion  planning  algortthu"'-wiil  have  to  find  k 
mutually  vertex  disjoint  paths,  one  connecting  "eaeh  s^  and  t^^, 
1  <  i  <  k.   QED 

(Remark:  k  must  grow  with  lv|). 
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Appendix  III 


Proof  of  (P7) 

We  shall  transform  any  instance  I  of  INDEPENDENT  SET  IN  CUBIC 
PLANAR  GRAPHS  [GJ,  pg.  194]  to  an  instance  F^  of  (Pf ) ,  such  that  the 
motion  in  F-j.  can  be  planijed  iff  the  graph  in  I  has  an  independent  set. 

Let  G  =  (V,E)  b^^any  instance  of  a  connected  cubic  planar  graph  (a 
small  modiflcatiRH  of j  ,  th^  construction  can  extend  to  non-connected 
graphs).  Let  k  ^e^an  integer,  <  |v|.  We  augment  G  as  follows:  Each 
vertex  v  of  G  is  connected  by  a  new  edge  to  a  new  vertex  u  .  The 
degree  of  each  u  is  one.  We  then  put  k  pebbles  in  k  of  the  new 
vertices.  Let  FP  =  V  and  let  us  assume  the  existence  of  a  motion 
coordination  algorithm  A  which  achieves  the  goal  of  the  problem.  We 
then  modify  A  so  that  it  moves  each  pebble  alone  from  its  initial  to 
its  final  position  (call  it  x)  and  then  it  moves  it  to  u  ,  Other 
pebbles  cannot  block  its  motion  as  soon  as  it  is  restricted  on  vertices 
of  G.  The  final  placement  of  vertices  clearly  defines  an  independent 
set  of  G,  of  size  k.  Conversely,  if  G  has  an  independent  set  S  of  size 
k,  then  there  is  an  easy  motion  planning  algorithm  which  moves  the 
pebbles  to  that  set:  It  moves  each  pebble  from  its  initial  position  to 
some  vertex  of  S  (by  connectivity  of  G  all  vertices  of  S  are  reachable) 
and  then  moves  the  pebble  to  u  so  that  v  can  be  used  for  motion  of 
other  pebbles  to  other  vertices  of  S.  Finally,  V  v£  S,  u  will  have  a 
pebble.  Then  the  algorithm  just  moves  each  pebble  form  u  to  v.  The 
time  complexity  for  planning  the  motion  given  that  S  is  known  is  O(nK). 


.A.^     A 


FIGURE   A.l 


FIGURE    A. 2 
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